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Abstract 



Electron tunneling through mesoscopic metallic grains can be treated per- 
turbatively only provided the tunnel junction conductances are sufficiently 
small. If it is not the case, fluctuations of the grain charge become strong. 
As a result (i) contributions of all ~ including high energy - charge states 
become important and (ii) excited charge states become broadened and es- 
sentially overlap. At the same time the grain charge remains discrete and the 
system conductance e-periodically depends on the gate charge. We develop a 
nonperturbative approach which accounts for all these features and calculate 
the temperature dependent conductance of the system in the strong tunneling 
regime at different values of the gate charge. 

PACS numbers: 73.40 Gk, 73.40 Rw 
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Coulomb effects may liave a strong impact on electron transport through small tunnel 
junctions and metallic grains Provided the resistance of tunnel junctions Rt between 

the grain and the lead electrodes is large Rq = 7r/2e^ ~ 6.5 Kf2 many features of single 

electron tunneling are well described within a simple perturbation theory in a dimensionless 
junction conductance at = Rq/Rt This theory uses the concept of discrete charge states 
of the grain and describe the system dynamics in terms of occupation probabilities of these 
charge states. 

In the limit ^ 1 nontrivial features appear only in the vicinity of the Coulomb blockade 
threshold in which case two charge states become nearly degenerate and the perturbation 
theory fails. Nonperturbatively this problem has been treated in Refs. and again the 

physical picture of electron tunneling via discrete charge has been exploited. 

If the conductance of tunnel junctions is not small at ^ 1 the problem turns out to 
be more complicated. Indeed, already making use of a simple perturbative formula for the 
inverse lifetime of the excited grain charge state Q = Qq > e/2 at T = (see e.g. [0,0) 
r = 2ate{Qo — e/2)/7iC one can immediately conclude that for > 1 broadening of the 
excited charge states Q ^ e due to strong quantum fluctuations of the charge is of the order 
of the spacing between them F ~ Ec- Thus charge levels overlap and the very concept of 
tunneling via discrete charge states with given energies becomes illdefined for such values of 
a. The finite T effect makes this overlap even more pronounced. 

In this Letter we propose a theoretical approach which allows to obtain a quantitative 
description of electron transport through mesoscopic metallic grains in the strong tunneling 
regime. We reformulate the problem in terms of the variable canonically conjugated to 
that of the charge, analyze its quantum dynamics and obtain an expression for the system 
conductance valid for all values of the gate charge and practically all experimentally relevant 
values of temperature. 

We will consider a standard model for a SET transistor: a small metallic grain is em- 
bedded between two bulk electrodes and connected with them via tunnel junctions with 
resistances Rl and Rr and capacitances Cl and Cr respectively for left and right junctions. 
A gate voltage Vg is applied to the grain via a capacitance Cg, a transport voltage between 
two electrodes is equal to V^. We also assume that the impedance of the external circuit is 
much less then quantum resistance. This system can be described by the Hamiltonian 

H={q-Qg)V2C + HL + HR + Hg + HT (1) 

where q is the operator of the grain charge, Qg = ClVl — CrVr + CgVg is the (noninteger) 
external charge, Vl,r. = VxRl,r/ {Rl + Rr) and C = Cl + Cr + Cg. The terms Hh = 
J2k ^khdkh^kh describe the kinetic energy of noninteracting electrons in the left {h = L) and 
right [h = R) electrodes and in the grain {h = g), whereas the term 

Ht= J^^hatgttk'heyipi-i'^h/'^) + c.c. (2) 

h=L,Rk,k' 

takes care about electron tunneling between the electrodes and the grain (the tunneling 
matrix elements T^',^ multiplied by the densities of states yield the junction resistances 
Rl,r = 4:716^ NL R{0)Ng{0)\Ti R\'^ ). The junctions phase operators can be expressed as 
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2eCR,z. 2eCgt ( Cr-Cl V. 

Cl + Cr " ^ Cl + Cr + C, \Cr + Cl2 ' 



where (p is the "phase" of the grain, or, more exactly, is the operator canonically conjugated 
to the grain charge q: [(^, = 2ei. 

After a standard procedure of averaging over the electronic degrees of freedom (see 
e.g. [0]) one can reformulate the problem in terms of the reduced density matrix p{ip,ip') 
which depends only on the phase variable (p. If the charge varies continuously everywhere 

in the system the density matrix pdf, f') is nonperiodic in (y9 and obeys a standard 

+00 

normalization condition / dippc{(p, (f) = 1. In our physical situation, however, the charge 

—00 

on the grain is quantized in units of the electron charge e. In this case the phase variable 
is compact (i.e. the states (p and + 47r are equivalent) and the density matrix obeys the 
conditions 

Pd(V3i + 47m, + 47rm) = exp (^'^^^^^{n - m)^ pdivi, (3) 



2-K 

dippd{ip,ip) = 1. (4) 
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Let us now introduce a nonperiodic in ip density matrix p{ipi,ip2)'- 

pdi'Pi, V2) = - m)] p{ipi - Ann, - 47rm), (5) 



which satisfies the following normalization condition 

j dip exp -71^ p[ip — Ann, (p) = 1 (6) 

^ -00 

The matrix p (|^) obeys the same equation of motion as the density matrix pc describing 
the continuous charge distribution in the system. If we assume that our system is ergodic 
we can immediately establish the connection between these density matrices. Indeed, at 
sufficiently large time the solution of the linear equation of motion for a dissipative ergodic 
system acquire a unique form irrespectively to a particular choice of the initial conditions. 
Thus in the long time limit the two solutions of this equation may differ only by a constant, 
which can be fixed with the aid of the normalization condition and we get ai t ^ 00 

~(. N Pc{t,ipi,ip2) , . 

p{t, y^i, ^2) = — — . (7) 

E / dcpexpU^^n) pdt,(p-Ann,ip) 

n -00 ^ ^ 

The equation (^ can now be used for evaluation of the expectation value of an arbitrary 
operator A{ip), which is 47r-periodic in ip. With the aid of (^,0) after a simple algebra we 
obtain 
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E (i(^) exp (^M%^n) 



A) = / d^A{^)pa{v,^) = ^^—- x\ (8) 



This is one of the main results of the present paper. It estabhshes a straightforward con- 
nection between the expectation values for an operator of any physical quantity calculated 
for discrete and continuous charge distributions. 

Let us now turn to a calculation of the tunneling current through a SET transistor. We 
first obtain a formal expression for the expectation value of the current operator in our 
system and then evaluate it with the aid of the equation @. The first part of this program 
will be carried out within the quasiclassical Langevin Equation approach [^,0 developed 
under the assumption that fluctuations of the phase variable are (in some sense) small. This 
is a suitable assumption as long as fluctuations of the charge are large. Expressing the kernel 
of the evolution operator in terms of the path integral on a real time Keldysh contour and 
calculating this integral within the quasiclassical approximation (see for further details) 
we obtain 



(9) 



where ipL^R/2e and qi^R define respectively fluctuating voltages and currents across the left 
and the right junctions, iL,Ri,2 are Gaussian stochastic variables describing the shot noise 
in these junctions and obeying the conditions 
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/ 27r * 


k2t) 



As the external impedance is negligible, the phases (Pl,r are linked to the transport voltage 
Vx by means of an obvious equation ipi + ^r = SeV^. Further relations between the phase 
and the charge variables are defined by the charge conservation law Ql — Qr = Cgipg/2e 
and the Kirghoff's law (phj^e — 'pR/2e = 2Vg — 0g/e, (pg is the gate capacitance phase 
defined analogously to (Pl,r- Combining all these equations with @ after averaging over 
the stochastic variables ^ we arrive at the expression for the current in our system: 



, , K RlUl)+RrUr.^ ^ ^ 

To evaluate the average values in (^) we shall use the result (|]). Assuming that fluctu- 
ations of the charge are Gaussian the contribution of the n-th term to the expectation value 
(|) can be roughly estimated as: 

Ai^) exp (^'-^n) ) ^ exp (-^^^n^) . (11) 

Thus provided the charge fluctuations are not small ^ it is sufficient to leave only 
the terms with n, m = 0, ±1 in the expression (§). In this approximation we obtain 
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^R,l 



1 + 2 (cos 



(12) 



where (...) denotes the average with the density matrix pc describing the continuous charge 
distribution. Making use of the equations (ll^), dH) and assuming the phase fluctuations to 
be small ^ vr in the limit of small transport voltages we arrive at the expression for the 
linear conductance 



{Rl + Rr)G{T) = 1 - /(T) - (7(T)e-^(^) cos 



2nQ, 



(13) 



where Qg = CgVg. We define at = 7r/2e^-Ro, V-Ro = 1/-Rl + 1/-Rij and 



f{T) 



I 

2^t 



2mK 



+ ^ 1 + 



2atEc 

2^^ 



TC 



(14) 



2TT^{5q\T)) _ TT T xcoth( 



^(^) . . __ J ^ (15) 



-oo 



^/ X TT^T , / {2atEc\ , / 2atEc 



^?(^) = V / (^h^Lrt) ^ (^^^^ (cosh«(t) - 1) + ^ir(t)sinhn(t)) . (16) 



Here '^{x) is the logarithm of the gamma-function, 7 = 0.577... is the Euler constant and 

K{t) = Roe{t){l - exp(-t/i?oC)), (17) 
m D f A coth^sina;t 

Note that the value i^(0) oc (5g^(0)) ( plSf ) diverges logarithmically at high frequencies. As 
in the small voltage limit the Langevin equation approach does not work at very low T 
(see below), in order to define -F(O) (or, equivalently, the high frequency cutoff for ([T5|)) we 
should make use of a more rigorous technique. For = we find 

/ ^^27rg\ _ j dippcc^{AT[ + ip,(p) 2'KQg 

Peq{<^, <^') is the equilibrium density matrix of our system. In the limit > 1 this matrix 
was evaluated by means of various nonperturbative approaches In the leading order 

in e~^"* all these approaches yield /'^'^^°q(^''+''^''^) _ g-2at_ xhus we obtain 

F{0) - 2at. (18) 
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This equation completes our results. In the limit of large T the Q^-dependent part of the 
conductance vanishes and we find the asymptotic behavior: 



(fi„ + fl.)G(r) = l-§ + 5^a,(f )'-... (19) 

At lower temperatures the conductance suppression due to charging effects becomes more 
pronounced (fig. 1). Furthermore, by changing the gate charge Qg it becomes possible to 
e-periodically tune the value of G. The minimum and maximum conductance values (fig. 
1) correspond to = and Qg = e/2 where the Coulomb barrier for electron tunneling 
reaches respectively its maximum and minimum values. The modulation of G with Qg also 
becomes more pronounced as the temperature is lowered (fig. 2). 

Finally, let us formulate the validity condition for our results. Analogously to ||^ we find 
that the phase fluctuations are sufficiently small provided T ^ Ec for <^ 1 and 

T :^ atEc exp{-2at) (20) 

for at ^ 1. Another our assumption {Sq'^) ^ is also well justified for such values of T. 
According to ( pO] ) the validity domain of our analysis expands rapidly with increasing at- 



E.g. for the parameters of figs. 1,2 the condition (pO|) yields T ^ 20 mK. And indeed our 
results show a very good agreement with a preliminary experimental data of the Saclay 
group [1^ down to such small values of T. 
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FIGURES 

FIG. 1. Maximum and minimum conductance versus temperature {Ec = 0.715K, at = 2.12) 
FIG. 2. Dependence of conductance on gate voltage {Ec = 0.715K, at = 2.12) 
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